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Abstract 

This paper presents a new method for comput- 
ing acoustic signals from helicopter rotors in for- 
ward flight The aerodynamic and acoustic 
solutions in the near field are computed with a 
finite-difference solver for the Euler equations. A 
nonrotating cylindrical Kirchhoff surface is then 
placed around the entire rotor system. This Kirch- 
hoff surface moves subsonically with the rotor in 
forward flight The finite -difference solution is 
interpolated onto this cylindrical surface at each 
time step and a Kirchhoff integration is used to 
carry the acoustic signal to the far field. Com- 
puted values for high-speed impulsive noise show 
excellent agreement with model-rotor and flight- 
test experimental data. Results from the new 
method offer high accuracy with reasonable com- 
puter resource requirements. 

Introduction 

In addition to the desire for high aerodynamic 
performance, modern helicopter designs also aim 
for low rotor noise. This is particularly important 
for civilian helicopters that operate near heavily 
populated areas. 

There are two main types of noise that cause 
problems for helicopters. The first type is noise 
that is due to the interaction of the rotor blades 
with their vortical wake systems. This type of 
noise is called blade-vortex interaction, or BVT, 
noise. The second type of noise is called high- 
speed impulsive, or HSI, noise. It is characterized 
by a strong acoustic disturbance that occurs over 
a very short period of time. Impulsive noise is gen- 
erally associated with high tip speeds and advanc- 
ing tip Mach numbers greater that 0.9. 

Accurate prediction of rotor noise is essential 
for its control. The most commonly used noise pre- 
diction techniques are based on the Ffowcs Will- 
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iams and Hawkings equation [1]. This approach 
contains terms that model three different compo- 
nents of rotor noise. The first two components are 
thickness noise and loading noise. These are com- 
puted from the linear superposition of integrated 
monopole and dipole sources over the surface of 
the blade. The third term is a nonlinear quadru- 
ple integral that is much more difficult to evalu- 
ate and typically neglected. Examples of this type 
of acoustics model are given in Refs. [2-3]. 

The difficulty in modeling the nonlinear qua- 
drupole term is the main drawback with acoustics 
models that are based on the Ffowcs Williams and 
Hawkings equation. Without this term, the acous- 
tic signals in the far field are typically underpre- 
dicted as shown in Ref. [4]. 

Improved accuracy has been obtained with 
alternate methods that are based on nonlinear 
computational fluid dynamics (CFD). For exam- 
ple, Baeder [4] solved the Euler equations to 
model the acoustics of rotor blades both in hover 
and in forward flight. Acoustic solutions were 
obtained at distances of up to 3.5 radii from the 
rotor hub. The problem with this approach is that 
the demand for computer resources increases 
exponentially as the solution domain is extended 
beyond the rotor blade. It is not currently practi- 
cal to propagate helicopter acoustic waves much 
beyond 3 rotor radii without excessive numerical 
dissipation. 

A third approach to rotor acoustic prediction 
uses the combination of a CFD method close to the 
rotor blade and a linear Kirchhoff integral for- 
mula to carry the acoustic solution to the far field. 
The Kirchhoff integral approach, such as that in 
Ref. [5], integrates a known pressure field over a 
prescribed surface and then propagates this sig- 
nal to arbitrary distances from the rotor blade. 
The CFD method accurately captures the tran- 
sonic flow nonlinearities close to the blade, while 
the Kirchhoff integral scheme is computationally 
more efficient over large distances. 

Two types of hybrid CFD/Kirchhoff methods 
have been demonstrated for rotary- wing applica- 
tions. The two methods differ as to whether or not 
the Kirchhoff surface rotates with the blade. 
Lyrintzis et al. [6-8] use a Kirchhoff surface that 
rotates with the blade. Pressure data on the 



Kirchhoff surface are computed from a numerical 
solution of the full -potential equations. The Kirch- 
hoff integral for the moving surface is computed 
with the formulation given in Ref. [5], This 
approach has the advantage that the Kirchhoff 
integral uses the same computational mesh as the 
CFD calculation. Also, the Kirchhoff surface can 
be positioned to minimize numerical dissipation in 
the CFD solution. 

However, this rota ting-surface Kirchhoff formu- 
lation has one major problem. The Kirchhoff inte- 
gral in Refs. [5-8] assumes that the Kirchhoff 
surface moves subsonically. Evaluation of the inte- 
gral for a surface that moves supersonically is 
much more difficult and has not been successfully 
applied to a helicopter problem. 

Restriction of the Kirchhoff integral to subsonic 
surface motion can be a major problem for high- 
speed rotary- wing applications. This is because the 
rotating velocity increases as £2 r , where £2 is the 
blade angular velocity and r is the distance in the 
plane of the rotor from the hub to a point on the 
Kirchhoff surface. For the high-speed test case 
computed later in this paper, a rotating Kirchhoff 
surface must be located less than 1.4 chords from 
the tip of the blade in order to ensure subsonic 
motion. This location may be too close for accurate 
acoustic predictions because the strong aerody- 
namic shock on the blade surface creates nonlin- 
earities near the tip. 

The acoustic prediction methods in Refs. [3,9] do 
not have this problem with supersonic motion of 
the Kirchhoff surface. This is because the integral 
evaluations take place on a nonrotating surface. A 
coordinate transformation is used to interpolate 
the near-field CFD solution onto a nonrotating, 
cylindrical Kirchhoff surface. The challenge with 
this method is whether the CFD solution can carry 
the acoustic signal out to the Kirchhoff cylinder 
with low numerical dissipation. However, there is 
no constraint on the radial location of the surface 
as long as it completely encloses the rotor blades. 

Baeder et al. [3] have used a structured-grid 
Euler CFD solver and the nonrotating Kirchhoff 
formulation to compute HSI noise for a hovering 
rotor. Strawn et al. [9] have used a similar scheme, 
but their Euler solver used unstructured, solution- 
adaptive grids to improve the resolution of the 
near-field acoustic signals. Both methods have 
shown excellent agreement with experimental 
data for HSI noise from hovering rotors. 

Both the Baeder et al. [3] and Strawn et al. [9] 
Kirchhoff integration schemes were limited to hov- 
ering rotors. Rotors in forward flight are signifi- 
cantly more difficult because the CFD solutions 
must be computed and stored at many time steps 


for the unsteady flowfield. The purpose of the cur- 
rent work is to develop a new CFD/Kirchhoff 
scheme for predicting helicopter acoustics in for- 
ward flight. The new formulation uses a nonrotat- 
ing Kirchhoff surface that moves subsonically with 
the rotor hub. The Kirchhoff integral is evaluated 
with the formulation given in Ref. [5]. Develop- 
ment of appropriate data structures and efficient 
interpolation schemes are the primary tasks in 
this effort 

This is the first time that a nonrotating Kirch- 
hoff surface has been used to compute the acous- 
tics from rotors in forward flight Results from the 
method are compared to experimental data for 
HSI noise. Solution accuracy is addressed and 
computed solutions exhibit minimal numerical dis- 
sipation. The overall computational efficiency of 
the method is also discussed. 

Near-Field CFD Solution 

The structured-grid Euler/Navier-Stokes solver 
called TURNS [10,11] is used to compute the aero- 
dynamic field close to the helicopter rotor. This 
CFD code solves the Navier-Stokes equations 
about rotating helicopter blades. Since viscous 
effects are minimal for the test cases considered in 
this paper, the TURNS code is run in an inviscid 
mode. 

The two computed cases in this paper were 
experimentally tested by Schmitz et al. [12]. They 
consist of a 1/7 scale research model of a US Army 
AH-1 helicopter with a blade aspect ratio of 9.22. 
Both cases have the same hover- tip Mach number 
of 0.665, with advance ratios of 0.258 and 0.348, 
respectively. 

Identical computational grids for the TURNS 
code have been constructed for both test cases. 
They consist of a series of 50 C-meshes that are 
stacked in the spanwise direction, with 20 located 
on the blade surface. Because the computed cases 
are nonlifting, the problem is symmetric about the 
plane of the rotor. This means that the solution 
need only be computed over half the computational 
domain. Each C-mesh is located along a constant 
radial line from the hub of the rotor and contains 
68 points in the chordwise direction with 48 points 
on the lower surface of each airfoil section. 35 
points are located in the direction normal to the 
blade surface. 

A view of the CFD mesh in the plane of the rotor 
is shown in Fig. 1. Note that the computational 
domain extends out to 2 rotor radii from the hub in 
the plane of the rotor. The outer boundary of the 
grid below the rotor blade is set at 1.5 radii. 
Between the blade tip and the outer spanwise 


boundary, the clustered region of the mesh is 
swept backwards In an effort to capture the acous- 
tic signal with minimal numerical dissipation. 

Similar solutions for these cases were computed 
on similar grids by Baeder [4], who also used the 
TURNS code. In that study, however, the acoustic 
field was computed directly with the CFD code and 
the Kirchhoff surface approach was not used. 
Because of this, Baeder’s CFD solutions covered a 
much larger computational domain in the span- 
wise direction than in the current work In spite of 
these differences, the acoustic solutions near the 
blade tip in this paper are virtually identical to 
those computed by Baeder. 

The TURNS code is first run in the quasi-steady 
mode to determine a starting solution at zero 
degrees azimuth angle. This requires approxi- 
mately 1800 iterations corresponding to 24 CPU 
minutes on a Cray C-90 computer. The unsteady 
time marching is then started with each time step 
corresponding to 0.25 degrees of blade azimuth 
angle. Approximately one hour of C-90 CPU time 
is required to complete a full 360° of rotor motion. 

The TURNS code has been modified so that the 
acoustic pressure, p , as well as its normal and 

temporal derivatives, p n , and p (} are computed at 

each time step. Note that the p t derivative must be 

computed in a nonrotating reference frame so that 
it is compatible with the nonrotating Kirchhoff 



Figure 1: View of the 3-D CFD grid in the plane of 
the rotor. 


surface. Pressure values that are located on the 
Kirchhoff surface are written out to a file at inter- 
vals of one degree of azimuthal angle. This time- 
dependent data base is later used to evaluate the 
Kirchhoff integral. 

Kirchhoff Surface Method 

It is not practical to continue the CFD solution 
to large distances in the spanwise direction. Large 
numbers of mesh points are required and the cal- 
culation rapidly becomes too large for existing 
computers. An alternate approach is to place a 
nonrotating cylindrical Kirchhoff surface around 
the rotor blades as shown in Fig. 2. Strictly speak- 
ing, the Kirchhoff surface should completely 
enclose the rotor blades, but the top and bottom 
surfaces are neglected for these computations. 
They are located so far above and below the rotor 
plane that their contributions to the far-field 
acoustics are typically very small. Most of the rotor 
noise is produced in the plane of the rotor. 

The Kirchhoff surface translates with the rotor 
hub when the helicopter is in forward flight. The 
acoustic pressure, p y at a fixed observer location, 

x, and observer time, /, can be evaluated by per- 
forming the following integration on the cylindri- 
cal surface: 



Figure 2: Schematic for the Kirchhoff surface 
integration. 
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This formulation is taken from Ref. [5]. It 
assumes that the Kirchhoff surface is moving with 

Mach number M . The distance between a point on 
the Kirchhoff surface and the observer is given by 
|r|, as shown in Fig. 2. Also note that the entire 
integral in Eq. (1) is evaluated at the time of emis- 
sion for the acoustic signal, x. 

The expressions for E j and E 2 are given as: 
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These expressions assume that the surface is mov- 
ing with steady translational motion. Additional 
terms that are required to account for unsteady or 
rotational motion are given in Ref. [5], The expres- 


sion for £ 2 in Eq* (3) has been modified from the orig- 
inal Farassat and Meyers' [5] formula by using the 
simplified expression found in Ref [13]. 

In the above equations, M n and M r are the com- 
ponents of M along h and r in Fig. 2. M t is the veloc- 
ity vector tangent to the Kirchhoff surface, and Vp 
is the gradient of the pressure on the Kirchhoff sur- 
face. The freestream speed of sound is assumed to be 
uniform at a^, and the angle, 0, is defined in Fig. 2. 

Evaluation of the integral in Eq. (1) at the emis- 
sion time requires a series of coordinate transforma- 
tions to properly access the CFD database on the 
Kirchhoff surface. These transforms can be 
described with the aid of Fig. 3. Fig. 3a shows the 
rotor blade and Kirchhoff surface at the time the 
sound reaches the observer. However, Eq. (1) 
requires that the pressures on the Kirchhoff surface 
be evaluated at the time they were emitted. At the 
time of emission, both the Kirchhoff surface and the 
rotor blade were in different locations. 

In order to find these locations, the delay between 
the observer time, f, and the emission time, x, must 
first be computed. This can be determined from Fig. 
3b by noting that the time it takes the acoustic sig- 
nal to travel from the Kirchhoff surface to the 
observer is equal to the time it takes the translating 
Kirchhoff surface to move the distance, d. These 



Figure 3: Two successive coordinate transformations are applied to evaluate the pressure data on the Kirch- 
hoff surface. 


times are given as |r| /a^, and , respec- 

tively. This leads to a quadratic Equation for the 
required time delay. One of the roots is nonphysi- 
cal and can be discarded. The locations of the 
Kirchhoff surface and the rotor blade at the emis- 
sion time can then be computed. 

Once the geometry is established at the emis- 
sion time, the acoustic pressures and their deriva- 
tives on the Kirchhoff surface are interpolated 
from the stored CFD database (see Fig. 3c). The 
CFD solutions are stored at discrete time steps on 
the Kirchhoff surface as two-dimensional quadri- 
lateral meshes. The required acoustic data are 
determined by linear interpolation in time and 
space. 

This procedure must be performed for every dis- 
crete integration point on the Kirchhoff surface. 
Topical grid sizes for this integral evaluation are 
1440 points in the azimuthal direction and 100 
points normal to the plane of the rotor. The azi- 
muthal points are equally spaced to enclose both 
rotor blades. The vertical extent of the Kirchhoff 
surface is ±1.5 rotor radii, the same as that for the 
CFD grid. The point spacings in this direction are 
exponentially stretched from the plane of the rotor 
to the outer boundaries. 

Results 

The new hybrid CFD/Kirchhoff method has 
been used to predict the acoustic signals from two 
model-rotor wind-tunnel cases described in Ref. 
[12]. These experiments recorded acoustic signals 
from a 1/7 scale model of the Army’s AH-1 helicop- 
ter main rotor. Microphones were placed at several 
fixed locations around the rotor system. Scaled 
acoustic data from flight tests is also reported in 
Ref. [12]. 

The two computed test cases both a have hover- 
tip Mach number, M (ip , equal to 0.665. The 

advance ratios, p, for the low-speed and high- 
speed cases are 0.258 and 0.348, respectively. The 
rotor thrust coefficient is the same for both cases 
and is equal to 0.0054. These rotor blades have 
symmetric airfoil sections with a thickness-to- 
chord ratio of 0.0971 and an aspect ratio of 9.22. 
The primary noise-generation mechanism in both 
cases is high-speed impulsive (HSI) noise. 

In spite of the fact that the model rotor experi- 
ments have a significant amount of thrust, the 
computations in this paper are for nonhfting rotors 
with the collective pitch set to zero. This simplifies 
the analysis since the rotor wake does not have to 
be modeled in the CFD solution. In general, this is 


a complicated task and is not the focus of this 
paper. 

The justification for neglecting the rotor thrust 
is that HSI pressure signals in the plane of the 
rotor are generally insensitive to thrust. This has 
been experimentally documented by Schmitz et al. 
[12,14]. The nonlifting assumption in the analysis 
has little effect on the computed results as long as 
acoustic comparisons are restricted to the plane of 
the rotor. Acoustic predictions that are out of the 
rotor tip path plane will require realistic rotor 
wake models in the CFD solutions. This is particu- 
larly true for cases with blade-vortex interactions. 

The Kirchhoff surface is located at 1.39 rotor 
radii for both computations. This location is far 
enough from the blade tip that nonlinear transonic 
effects are small, but close enough so that numeri- 
cal dissipation does not degrade the CFD solution. 
Both Baeder et al. [3] and Strawn et al. [9] have 
investigated the choice of Kirchhoff surface loca- 
tions for high-speed hovering rotors. Their results 
show grid independence for Kirchhoff surfaces at 
1.4 radii. 

If comparisons to experiment are restricted to 
the plane of the rotor, then the Kirchhoff integra- 
tion is symmetric about this plane. This means 
that the integral in Eq. (1) need only be computed 
over half of the Kirchhoff surface. The resulting 
pressure can then be doubled to account for the 
remainder of the integration. As such, the numeri- 
cal integration on the Kirchhoff surface consists of 
1440 equally-spaced points around the azimuth 
and 50 unequally-spaced points in the lower half- 
plane of the rotor blade. 

With this grid, the Kirchhoff integration in Eq. 
(1) requires about 30 CPU seconds on the Cray C- 
90 for each evaluation of the observer pressure. 
Most of this time is spent performing interpola- 
tions in the CFD database. There is a potential for 
significant speedup if these interpolations can be 
performed more efficiently. 

Typical acoustic pressure signals in this paper 
consist of approximately 24 evaluations of the 
observer pressure. This requires a total CPU time 
of 12 minutes to obtain the complete acoustic sig- 
nal at each observer location. At large distances 
from the rotor blade, this is orders of magnitude 
less than the time that would be required to com- 
pute a pure CFD solution for the same location. 

Figure 4 compares the computed and experi- 
mental results for the low-speed case. This case 
has an advancing-tip Mach number, M af} of 0.837 

and the computations show significant transonic 
flow at the blade tip. This transonic flow is limited 
to the blade and does not connect to the far-field 
region of supersonic flow relative to the blade. The 






Figure 5: Acoustic pressure comparisons for M = 0.896. 


angle, where the surface shock at the blade tip is 
the strongest. The Mach-one lines are drawn darker 
than the other contours. Note that the supersonic 
region on the blade surface almost connects with 
the supersonic region in the far field. When this 
phenomenon occurs, it is referred to as delocaliza- 
tion and the surface shock is free to propagate to 
the far field with very little dissipation. The acous- 
tic amplitude increases dramatically at the onset of 
delocalization. 

The delocalization phenomena is highly depen- 
dent on nonlinear transonic effects that occur near 
the blade tip. Fig. 6 shows that the flowfield is not 
quite delocalized, but it should be noted that this 
CFD solution was computed for a nonlifting rotor. If 
the cyclic pitch and wake effects were included in 
the computation, it is reasonable to assume that 
these might have some effect on the picture in Fig. 


6. If these effects caused the flowfield to delocalize, 
then the predicted peak negative pressures would 
be larger, and thereby show better agreement with 
the experimental results. 

Experimental data from the flight test are also 
shown for microphone 2 in Fig. 5. This flight-test 
data shows excellent agreement with the computer 
predictions but not with the model-scale results. 
The reasons for the discrepancies between model- 
scale and flight test are not known. Perhaps it is 
related to the sensitivity of the far-field acoustic 
pressures when the flowfield at the rotor tip is very 
close to delocalization. 

In spite of these differences, the computed 
results in Fig. 5 are significantly better than those 
shown by other researchers. Baeder [4] computed 
this case with a pure CFD approach and obtained a 
peak negative pressure of 270 Pa for microphone 


M > 1 



Figure 6: Mach contours at 105° for M a( = 0.896. 


location 2. It is likely that part of this underpredic- 
tion is caused by numerical dissipation in the 
numerical solution at 3.44 radii. The CFD mesh 
did not include the 6.88 radii location in that 
study. Baeder also mentions that linear methods 
based on the Ffowcs Williams and Hawkings equa- 
tion predict a peak negative pressure of only 210 
Pa for the same location. The current computation 
gives a peak negative pressure of 3 12 Pa compared 
to the experimental result of 397 Pa. 

Discussion 

The accuracy of the calculations can be 
addressed by examining both the fundamental 
approximations in the Kirchhoff formulation and 
the mesh independence of the computed results. 
The first assumption in the Kirchhoff* formulation 
is that the Kirchhoff surface moves through undis- 
turbed air. This is not entirely true for a helicopter 
because lifting rotor blades generate aerodynamic 
disturbances in their wake systems. The interac- 
tion of the acoustic signals with the rotor wake 
outside the Kirchhoff surface is not considered in 
the present method.. The effect of this approxima- 
tion should be small however. This is because a 
rotor in forward flight propagates most of its 
acoustic disturbances ahead of it while the wake 
system is left behind. Thus there is little interac- 
tion between the forward-radiating acoustic sig- 
nals and the rotor wake. 

Another basic assumption is that the speed of 
sound is constant outside the Kirchhoff surface. 


This means that the surface must be located far 
enough from the rotor blade to completely enclose 
all nonlinear transonic effects. We can check 
whether this condition has been met for the com- 
puted results shown in Figs. 4 and 5 by performing 
a new computation with a different Kirchhoff sur- 
face location. If the predicted acoustic pressures 
differ, then there may be some nonlinear effects 
that are not completely inside the surface. 

Figure 7 shows computed pressure signals for 
the high-speed test case at microphone location 2. 
The results differ in the positions of the Kirchhoff 
surfaces, s/R. The first prediction has a surface 
located at 1.39 radii and the second at 1.28. The 
two computed results are virtually identical, which 
indicates that the nonlinear effects are completely 
contained inside the Kirchhoff surface. 

Figure 7 also provides evidence of grid indepen- 
dence of the CFD solution. If numerical dissipation 
played a role in the CFD solution, then the far- 
field acoustic pressures would be affected by the 
location of the Kirchhoff surface. 

A final question involves the numerical resolu- 
tion of the Kirchhoff integral in Eq. (1). The calcu- 
lations in Figs. 4 and 5 used a 1440 x 50 mesh on 
the lower half of the Kirchhoff surface. This corre- 
sponds to a constant azimuthal resolution of 0.25° , 
and a minim um vertical resolution of 0.01 chords 
at the plane of the rotor. Figure 8 compares the 
original results for the high-speed case at micro- 
phone location 2 with those from a finer Kirchhoff 



Figure 7: Comparison of results for microphone 
location 2 using two different Kirchhoff surfaces. 




Figure 8: Comparison of results with two different 
mesh resolutions on the Kirchhoff surface. 

mesh. This finer mesh has a resolution of 2880 x 
80, where the minimum vertical spacing at the 
plane of the rotor has been decreased to 0.005 
chords. The two results are virtually identical, 
indicating that the computed values from the orig- 
inal Kirchhoff integration are mesh independent 

Summary 

This paper presents a new method for comput- 
ing far-field acoustics from helicopter rotor blades 
in forward flight A solution to the Euler equations 
accurately models the nonlinear effects near the 
blade surface and a Kirchhoff integration propa- 
gates the near-field acoustic signals to the far field 
in a computationally-efficient manner. The key to 
the Kirchhoff formulation is the use of a nonrotat- 
ing surface which ensures that its motion is 
always subsonic. 

Close to the blade tip, computed results with the 
new method compare favorably with experimental 
results and predictions from pure CFD methods. 
However, the major advantage of the new method 
over its pure CFD counterparts occurs for far-field 
calculations. CFD computations are generally lim- 
ited to computational domains of less than 3 or 4 
rotor radii. The combined CFD/Kirchhoff method 
can compute acoustic signals at arbitrary observer 
locations with minim al numerical dissipation. 

Although the method is demonstrated for cases 
with high-speed impulsive noise, it should also be 


applicable to cases with blade-vortex interactions. 

The primary challenge for such computations will 

be the accurate modeling of the rotor wake system 

in the CFD solver. 
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